Roughness of stylolites: a stress-induced instability with non local interactions 
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We study the roughness of stylohte surfaces (i.e. natural pressure-dissolution surfaces in sedimen- 
tary rocks) from profiler measurements at laboratory scales. The roughness is shown to be nicely 
described by a self-afBne scaling invariance. At large scales, the roughness exponent is ("1 ~ 0.5 
and very different from that at small scales where ^2 ~ 1.1. A cross-over length scale at around 
Ac — 1mm is well characterized and interpreted as a possible fossil stress measurement if related to 
the Asaro-Tiller-Grinfeld stress-induced instability. Measurements are consistent with a Langevin 
equation that describes the growth of stylolite surfaces in a quenched disordered material with long 
range elastic correlations. 

PACS numbers: 83.80.Ab, 62.20.Mk, 81.40.Np 



Stylolites are geological patterns that are very com- 
mon in polished marbles, a material largely used for 
floors and walls of buildings and monuments. They ex- 
ist in many sedimentary rocks such as limestones, sand- 
stones or evaporites They are rock-rock interfaces, 
that are formed during diagenesis and result from com- 
bined stress-induced dissolution and precipitation pro- 
cesses They exist on a very large range of scales, 
from micro- meters to tens of meters. They are often ob- 
served as thin irregular interfaces that look like printed 
lines on rock cuts which explain their name. At larger 
scales, they are planar structures that are typically per- 
pendicular to the tectonic load (i.e. lithostatic pressure). 
They might form complex network with or without con- 
nections. 

However, despite their abundance, stylolites are as 
mentioned by Gal et al. 5] "among the least well- 
explained of all pressure-solution phenomena". First they 
are complex 3D structures that are often only described 
from 2D cross-sections. Second, they develop in various 
petrological and tectonic contexts which lead to very dif- 
ferent geometries. Third they are often transformed and 
deformed because of post-processes like diagenesis and 
metamorphism that develop after their initiation. 

In this letter we propose a new mechanism for stylo- 
lites growth that is consistent with recent 3D roughness 
measurements of a stylolite interface. The first part of 
the letter deals with the topography measurement and 
its description in terms of scaling invariance, namely self- 
affinity. The second part is devoted to a physical model 
of the stylolites roughening. It is based on a Langevin 
equation that accounts for a stress-induced instability in 
a quenched disorder with capillary effects and long range 
elastic interactions. 

The roughness measurement has been performed on a 
stylolite interface included in a limestone sample from 
Juras Mountains in France (see Fig. The sample has 
been collected in a newly open quarry, thus preserved 




FIG. 1: Picture of the stylolite surface that has developed in a 
limestone from Juras Mountains (France). Magnitude of the 
peaks are typically of the order of 2 mm. The width of the 
sample is 4 cm. 



from late breakage and chemical erosion. The opening 
procedure was possible for this sample because of the ac- 
cumulation of undissolved minerals like clays that form 
a weak layer along the stylolite interface. The concen- 
tration of these minerals provides an estimate of the cu- 
mulative strain that the sample has undergone 0]. As 
shown in Fig. ^ peaks along the interface are randomly 
distributed in space and of various sizes (up to several 
millimeters). Locally slopes along the topography might 
be very important which makes the roughness measure- 
ment difficult. 

We used two different profilers to sample the stylolite 
roughness. First, with a mechanical profiler 0,13 we ex- 
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FIG. 2: A ID profile obtained by a mechanical profiler (1030 
data points - Ax — 30^m). 



tracted four profiles of 1030 points each with a horizontal 
step of Ax — 30 fim. The specificity of the mechanical 
profiler is to measure the surface height from the contact 
of a needle onto the surface. The radius of the needle tip 
is 25 fim. The vertical resolution is 3 fim over the avail- 
able range of 5 cm. One of this profile is shown in Fig. |21 
We compare the mechanical measurement to an optical 
profiling 1^. This technique is based on a laser triangu- 
lation of the surface without any contact with the sur- 
face. The laser beam is 30 /^m wide. Horizontal steps be- 
tween measurement points were Ax = Ay = 50 iim. The 
main advantage of this technique is the acquisition speed 
that can be significantly larger compared to the mechan- 
ical profiler since there is no vertical move. However, a 
successful comparison with mechanical measurement is 
necessary to ensure that optical fluctuations are height 
fluctuations and not material property fluctuations. 

Fig. O shows a gray level map of the stylolite surface 
heights as a high resolution image (600 x 600 pixels) that 
covers over a surface area of 3 x 3 cm^ . The vertical reso- 
lution is 2 fim over the 6 mm vertical range of the setup. 




FIG. 3: 2D map obtained by an optical profiler (600x600 
data points - Aa; = Ay — 50/im). Heights are gray encoded: 
white corresponds to the minimum height and black to the 
maximum height. 



We analyzed the height distribution in terms of self- 
affinity [10| which states that the surface remains sta- 
tistically unchanged for the transform: Ax X Ax, 
Ay — !■ A Ay, Az Az, where A can take any 

real value. The exponent C is the so-called roughness 
exponent. A ID Fourier spectrum of a self-afflne pro- 
file is shown [l(l | to behave as a power law with a slope 
— 1 — 2C, and provides an estimate of the roughness expo- 
nent ^. We applied this technique whose reliability has 
been tested [ll|, to profiles extracted from the optical 
map shown in Fig. 13 either along x or y directions. Re- 
sults are plotted in Fig. 01 The figure shows with a thick 
solid line the y-direction average spectrum of 600 profiles 
extracted along the x-direction. The spectrum clearly 
exhibits two regimes. At small wavenumbers (i.e. large 
length scales), a power law behavior is observed with a 
slope in the log-log plot of -2 which is consistent with a 
roughness exponent of Ci = 0.5. At large wavenumber 
(i.e. small length scales), a second power law behavior 
is observed with a larger slope (-3.2) in agreement with 
a roughness exponent C2 = 1.1. The crossover length 
scale is sharp and defines a characteristic length scale 
Ac ~ 1mm. The average spectrum of profiles taken along 
the perpendicular direction (y-direction) provides a very 
consistent result and confirms the isotropy of the surface 
in terms of scaling invariance. We checked in Ref. @ 
that mechanical profiles show the same properties spe- 
cially at large wavenumbers since they are sampled at a 
higher resolution and allow an extended description of 
the unusual high roughness exponent C2 = l-l- We also 
checked that another analysis techniq ue, namely the Av- 
erage Wavelet Coefficient technique [IjI, was providing 
very consistent results. We also show in Ref. |^ that 
this behavior can be observed in several other stylolite 
samples. The difference between samples comes from 
the crossover length scale that can be larger or smaller, 
changing drastically the aspect of the surface. However, 
the two power law behaviors are conserved with the same 
values of the roughness exponents. 

The second part of the letter is devoted to a model- 
ing of the stylolite roughening in order to understand 
the origin of the self-affine behaviors and the parame- 
ter sensitivity of the characteristic length Ac. As usually 
done, we reduce the problem to a 2D situation z{x), i.e. 
invariance along the y-direction. 

We propose a phenomenological a ppr oach that follows 
the work reported in Kassner et al. |l3 who studied the 
dynamics of a strained solid in contact with its melt and 
more specifically the Asaro-Tiller-Grinfeld instability. In 
the case of a solid/fluid interface, the chemical potential 
between the solid and the fluid at the boundary can be 
written as |5llT3|: 

Afi = n{Ue+-fK) (1) 

where Ue is the elastic energy per unit volume in the 
solid, 7 is the surface energy, k the curvature and fl a 
molecular volume. We have assumed that gravity effects 
are negligible. If a bulk diffusion holds in the fluid, the 
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FIG. 4: Averaged power spectrum of the topographic profiles 
extracted from the optical map of the stylolite surface. 



evolution of the interface is directly related to the chemi- 
cal potential: Vn = mA/i where w„ is the normal velocity 
and m is the mobility. 

Following the work of Grinfeld and Hazzledine [l5j |. 
we propose to extend the approach to the case of a 
solid/solid interface. Moreover, we shall also assume the 
presence of a fluid in the pores. Porosity of the mate- 
rial is supposed to be sufficiently high for having both a 
bulk diffusion in the fluid and capillary effects. Doing so, 
Eq.(^ holds. Since a complete match between both solid 
surfaces is assumed, the normal velocity of the interface 
is now driven by the relative dissolution between the two 
solids. 

The difficult task is to expand the local elastic energy 
along the interface that takes into account the surface 
corrugations. Using the representation theorem for the 
displacement field and the Hook's law, it is possible to 
include the influence of non local elastic interactions in 
the limit of small perturbations of the interface through 
the elastic Green function as 0, 0| : 



a{x) = (To 



z{x') — z{x) 
{x' — x)^ 



dx' 



(2) 



where (t{x) is the local stress along the interface between 
the solid, co is the average uniaxial external stress or an 
effective stress that includes the fluid pressure, and PV 
is the principal value. According to Ref. the elastic 
energy can be approximated for plain strain by: 
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where E is an effective Young's modulus and v an effec- 
tive Poisson coefficient. 

The last aspect of the modeling concerns a descrip- 
tion of the noise that exists in the chemical potential be- 
cause of material fluctuations. We assume that this noise 



fluctuates on scales significantly smaller than the scales 
where elasticity and capillary effects are considered as for 
a Langevin approach However, unlike classical use 
of a Langevin equation, we assume here that the noise 
r]{x,z(x)) is quenched to represent spatial material fluc- 
tuations, and known to have a strong influence on the 
scaling properties. For instance, for the Edwards and 
Wilkinson (EW) problem 18, 19], the roughness expo- 
nent is 0.5 in the case of an annealed noise and 1.2 for a 
quenched noise. 

In the present framework, the interface evolution is 
described by the following equation: 
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which includes three driving forces: the long range elastic 
interactions, the local capillary effect which has a stabi- 
lizing influence, and the noise fluctuations. A crossover 
length scale is obtained by the balance of the average 
elastic and capillary terms in Eq.l^J l^j : 



A* 



jE 



(5) 



Typical values for limestones are E = 8 ■ 10^° N/m^, 
ly = 0.25, and 7 = 20 J/m^. If we consider a stress at 
2 km depth (cto = 40 MPa), we obtain A* w 1 mm. We 
propose that the A* length scale is the crossover length 
scale Ac observed in Fig. ^ 

Actually, the mechanical regime (A ^ Ac) and the cap- 
illary regime (A <C Ac) have both been extensively ex- 
plored but independently. Indeed, the capillary regime is 
nothing else than the EW problem in a quenched noise. 
Roux and Hansen ^19] have shown in this case that the 
interface is self-affine with an exponent (^2 ~ 1-2. In the 
mechanical regime, Eq.(@l) is reduced to the quasistatic 
propagation of an elastic line in a disordered material. 
This has been largely studied [lO, |^ 112 and the rough- 
ness exponent has been shown to be Ci « 0.4. 

We performed a numerical simulation of Eq.l^J us- 
ing an even driven algorithm [23, where there is no 
time evolution. The periodic interface is discretized in 
N=2048 cells. At start, the interface is flat. At each 
step, we are searching for the cell that exhibits the max- 
imum speed dz/dt according to Eq.Q). This cell is then 
advanced by a random amount dz uniformly sampled 
in [0,1]. The local fluctuation of the chemical poten- 
tial ri{x^z{x)) is updated from a prescribed distribution 
chosen here to be uniform between [0, 1]. Doing so, we 
are always dealing with the most unstable cell and let the 
interface grow there. 

After a transient regime, we observed a stationary 
width of the rough interface that develops. The power 
spectrum of the interface position P{k), averaged over 
500 samples, is shown in Fig.[Sl It reveals a behavior very 
similar to the stylolite measurements: A characteristic 
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FIG. 5: Averaged power spectrum of the topographic profiles 
obtained from the 2D styloUte modeUng. 



length scale defined as the cross-over between two power 
law regimes. The crossover length scale A* is controlled 
by the balance between the magnitude of the mechanical 
and capillary effects and is compared to the charateristic 
length scale A,;. As expected, the mechanical regime with 
an exponent (^i « 0.4, is dominating at small wavemim- 
bers {i.e. large length scales). On the contrary, at large 
wave numbers, the capillary regime dominates with a 
roughness exponent close to C2 = 1-2. Roughness expo- 
nents obtained from the modeling are slightly different 
from the measurements. This might be explained by the 
dimension difference: experimental surfaces are full 3D 



interfaces, on the contrary the model is 2D. A complete 
3D numerical modeling is on-going to account for this 
situation. However, 3D computations are much heavier 
because of the long delay to reach the stationary regime 
of the interface growth. 

In conclusion, we presented a quantitative description 
of stylolite interfaces that is consistent with a model of 
interface growth. The experimental measurement is a 
high resolution profiling of a 3D stylolite topography. We 
show that the surface is self-afJine but with two regimes. 
At small scales, the roughness exponent is unexpectedly 
high, (2 — 1.1, and consistent with a capillary dominated 
regime. At large scales, the stylolites morphology is con- 
trolled by mechanical effects, i.e. long range elastic stress 
redistribution. In this case the roughening is important 
with a low roughness exponent Ci = 0.5. The modeling 
is based on the description of a stress-induced instability 
previously reported as the Asaro-Tiller-Grinfeld instabil- 
ity. Such framework provides a prediction of a charac- 
teristic length Ac which is the crossover between the two 
scaling regimes. It is important for geological implica- 
tions to note that the characteristic length Ac is very 
sensitive to the average stress ctq. Indeed, a measure- 
ment of Ac from today roughness profiling could provide 
an estimate of the stress magnitude during the stylolite 
growth, that is, in the past. Thus stylolites could be 
considered as stress fossils. 
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